Event Excess in the MiniBooNE Search for v u — > v e Oscillations 
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The MiniBooNE experiment at Fermilab reports results from a search for t> M — > v e oscillations, 
using a data sample corresponding to 5.66 x 10 20 protons on target. An excess of 20.9 ± 14.0 events 
is observed in the energy range 475 < E„ E < 1250 MeV, which, when constrained by the observed 
v^ events, has a probability for consistency with the background-only hypothesis of 0.5%. On the 
other hand, fitting for v^ — > v e oscillations, the best-fit point has a ^-probability of 8.7%. The 
data are consistent with i> M — > i> e oscillations in the 0.1 to 1.0 eV 2 Am 2 range and with the evidence 
for antineutrino oscillations from the Liquid Scintillator Neutrino Detector at Los Alamos National 
Laboratory. 
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The MiniBooNE experiment has published searches for reported in |14| . A search in antineutrino mode provides 



and v n 



oscillations, motivated by the 



LSND 3.8(7 excess of v e candidate events [1(. In the 
Vy, — > v e study, MiniBooNE found no evidence for an 
excess of v e candidate events above 475 MeV; however, 
a 3.0cr excess of electron-like events was observed be- 
low 475 MeV_J2|, |3j. The source of the excess remains 
unexplained |3j, although several hypotheses have been 
put forward [4l-ll3|. including, for example, anomaly- 
mediated neutrino-photon coupling, sterile neutrino de- 
cay, and sterile neutrino oscillations with CP or CPT 



violation. Initial results from the v., 
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a more direct test of the LSND signal, which was ob- 
served with antineutrinos. Due to limited statistics, the 
initial MiniBooNE D^ — > D e search was inconclusive with 
respect to two-neutrino oscillations at the LSND mass 
scale, although a joint analysis reported compatibility 
between the LSND, KARMEN [TUlU, and MiniBooNE 
antineutrino experiments [6[. In this paper, we report an 
updated analysis of the v^ —¥ v e search with 1.7 times 



more protons on target (POT) than reported in [14|. 
This analysis uses the same technique that was re- 



study were ported earlier [14| and assumes only v^ — > u e oscillations 



with no significant z/„ disappearance and no v„ oscilla- 
tions. In addition, no contribution from the observed 
neutrino mode low energy excess has been accounted for 
in the antineutrino prediction. These simplifications may 
change the fitted v^ — > v e oscillation probability by a to- 
tal of - 10%. 

The antineutrino flux is produced by 8 GeV protons 



from the Fermilab Booster interacting on a beryllium tar- 
get inside a magnetic focusing horn. Negatively charged 
mesons produced in p-Be interactions are focused in the 
forward direction and subsequently decay primarily into 
v^. The flux for neutrinos and antineutrinos of all fla- 
vors is calculated with a simulation program using ex- 
ternal measurements [171 ]. In antineutrino mode, the v^,, 
v e , and v e flux contaminations at the detector are 15.7%, 
0.4%, and 0.2%, respectively. The v^ flux peaks at 400 
MeV and has a mean energy of 600 McV. 

The MiniBooNE detector has been described in detail 
elsewhere [18|. The detector location was chosen to sat- 
isfy L[m]/E[McV] ~ 1, similar to that of LSND, which 
maximizes the sensitivity to oscillations at Am 2 ~ 1 eV 2 . 
The detector consists of a 40-foot diameter sphere filled 
with pure mineral oil (~CH2). Neutrino interactions 
in the detector produce final-state electrons or muons, 
which produce scintillation and Cherenkov light detected 
by the 1520 8-inch photomultiplier tubes (PMTs) that 
line the interior of the detector. 

The signature of v^ — > T> e oscillations is an excess of 
£/ e -induced charged-current quasi-elastic (CCQE) events. 
Reconstruction [19( and selection requirements of these 
events are identical to those of the previous neutrino and 
antineutrino mode analyses |3j, |l4| . The detector cannot 
distinguish between neutrino and antineutrino interac- 
tions on an event- by- event basis. To help constrain the 
v e /v e candidate events, a v^jv^ sample is formed. The 
separation of v^ from v^ in this large CCQE sample is 
accomplished by fitting the observed angular distribu- 
tion of the outgoing muons to a linear combination of 
the differing CCQE angular distributions for final state 
/i + and \i~ . Relative to the Monte Carlo prediction, the 
fi + yield required an increase of 1.20 to the rate of ir~ 
decays {v^), while the [i~ yield is 0.99 of its predicted 
rate. Overall, the normalization required a 13% increase, 
which is compatible with the combined neutrino flux and 
cross section uncertainties [20J|. A sample of 24,771 data 
events pass the v^ CCQE selection requirements. The 
neutrino and antineutrino content of the sample are 22% 
and 78%, respectively. 

The oscillation parameters are extracted from a com- 
bined fit to the v e /v e CCQE and v^/v^ CCQE event 
distributions. Any possible v^ — > v e signal, as well as 
some v e backgrounds interact through a similar process 
as Vp CCQE events and are additionally related to the 
T>^ CCQE events through the same 7r~ decay chain at 
production. These correlations enter through the off- 
diagonal elements of the covariance matrix used in the 
X 2 calculation, relating the contents of the bins of the 
v e CCQE and v^ CCQE distribution. This procedure 
maximizes the sensitivity to v^ — > v e _ oscillations when 



systematic uncertainties are included [21 



The number of predicted v e CCQE background events 
for different ranges of reconstructed neutrino energy 
(E® E ) is shown in Table [TJ The background estimates 



TABLE I: The expected (unconstrained) number of events for 
different E„ E ranges from all of the backgrounds in the v e 
appearance analysis and for the LSND expectation (0.26% os- 
cillation probability averaged over neutrino energy) ofv^ — > v c 
oscillations, for 5.66 x 10 20 POT. 



Process 



200 - 475 MeV 475 - 1250 MeV 



v„ & v„ CCQE 


4.3 


2.0 


NC 7T° 


41.6 


12.6 


NC A ->■ iV 7 


12.4 


3.4 


External Events 


6.2 


2.6 


Other i/ft & v^ 


7.1 


4.2 


v e & v e from jti* Decay 


13.5 


31.4 


Ve & v e from K Decay 


8.2 


18.6 


Ve & v e from K ( [ Decay 


5.1 


21.2 


Other v e & v e 


1.3 


2.1 



Total Background 

0.26% V„ -> Ve 



99.5 
9.1 



98.1 
29.1 



include both antineutrino and neutrino events, the latter 
representing 44% of the total background. The predicted 
backgrounds to the v e CCQE sample arc constrained by 
measurements at MiniBooNE and include neutral cur- 

(NC) 



rent (NC) 7r u events [22j, A — > Nj radiative decays, 
and external events from neutrino interactions outside 
the detector. Other backgrounds from mis-identified Vy, 
or i?y j23l |24| and from intrinsic v e and v e events from the 
7r — > fi decay chain receive the v^ CCQE normalization 
correction according to their parentage at production (7r + 
or n~). 

If the low-energy excess observed during neutrino- 
mode running [2j were scaled by the total neutrino flux, 
the expected excess for antineutrino mode running would 
be 12 events for 200 < E® E < 475 MeV. (These events 
are not included in Tables HI and [TTI and in Fig. [1] because 
they occur below 475 MeV and the origin of these events 
is unexplained.) 

Systematic uncertainties are determined by consider- 
ing the predicted effects on the v^ and v e CCQE rate from 
variations of actual parameters. These include uncer- 
tainties in the neutrino and antineutrino flux estimates, 
uncertainties in neutrino cross sections, most of which 
are determined by in situ cross-section measurements at 
MiniBooNE, and uncertainties in detector modeling and 
reconstruction. By considering the variation from each 
source of systematic uncertainty on the v e CCQE signal, 
background, and v^ CCQE prediction as a function of 
E® E , a covariance matrix in bins of E® E is constructed. 
This matrix includes correlations between v e CCQE (sig- 
nal and background) and v^ CCQE and is used in the x 2 
calculation of the oscillation fit. 

Fig. Q] (top) shows the E® E distribution for v e 
CCQE observed data and background. A total of 277 
events pass the v e event selection requirements with 
200 < E$ E < 3000 MeV, compared to an expectation of 
233.8±15.3±16.5 events, where the uncertainty includes 
both statistical and systematic errors, respectively. This 
corresponds to an excess of 43.2 ± 22.5 events. (In the 
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FIG. 1: Top: The E® E distribution for u e CCQE data (points 
with statistical errors) and background (histogram with sys- 
tematic errors). Bottom: The event excess as a function of 
E„ . Also shown are the expectations from the best oscilla- 
tion fit with E$ E > 475 MeV, (Am 2 , sin 2 20) = (0.064 eV 2 , 
0.96), where the fit is extrapolated below 475 MeV, and from 
two other oscillation parameter sets in the allowed region. No 
correction has been made for the low-energy excess of events 
seen in neutrino mode below 475 MeV. All known systematic 
errors are included in the systematic error estimate. 

previous neutrino run analysis, event totals were consid- 
ered in two energy regions: 200 - 475 MeV and 475 - 
3000 MeV, where the latter region was the energy range 
for the neutrino oscillation search. For the antineutrino 
data, the excess for 475 < E$ E < 3000 MeV is 24.7±18.0 
events.) In the energy range from 475 < E® E < 1250 
MeV, the observed v e events, when constrained by the 
D,j, data events, have a x 2 /DF = 18.5/6 and a proba- 
bility of 0.5% for a background-only hypothesis. (This 
compares to the 40% probability that is observed in neu- 
trino mode [3| for the same energy range.) DF is the 
effective number of degrees of freedom from frequentist 
studies. The number of data, fitted background, and ex- 
cess events for different E® E ranges are summarized in 

TABLE II: The number of data, fitted (constrained) back- 
ground, and excess events in the v e analysis for different Ej 
ranges. The uncertainties include both statistical and con- 
strained systematic errors. All known systematic errors are 
included in the systematic error estimate. 



E® B Range 


Data 


Background 


Excess 


200 - 475 MeV 


119 


100.5 ±10.0 ±10.2 


18.5 ±14.3 


475 - 675 MeV 


64 


38.3 ±6.2 ±3.7 


25.7 ±7.2 


475 - 1250 MeV 


120 


99.1 ±10.0 ±9.8 


20.9 ± 14.0 


475 - 3000 MeV 


158 


133.3 ±11.5 ±13.8 


24.7 ±18.0 


200 - 3000 MeV 


277 


233.8 ±15.3 ±16.5 


43.2 ±22.5 
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FIG. 2: The E V i S (top panel) and cos(t9) (bottom panel) dis- 
tributions for data (points with statistical errors) and back- 
grounds (histogram with systematic errors) for E„ B > 200 
MeV. 

Table IE 

Fig. [2] shows the observed and predicted event distri- 
butions as functions of reconstructed E vis and cos(t9) for 
200 < E$ E < 3000 MeV. E vis is the measured visible en- 
ergy, while 9 is the scattering angle of the reconstructed 
electron with respect to the incident neutrino direction. 
The background-only \ 2 values for the v e and v^ data 
are \ 2 /DF = 23.8/13 and x 2 /DF = 13.6/11 for E vis 
and cos(6), respectively. 

Many checks have been performed on the data to en- 
sure that the backgrounds are estimated correctly. Beam 
and detector stability checks show that the neutrino 
event rate is stable to < 2% and that the detector energy 
response is stable to < 1%. In addition, the fractions of 
neutrino and antineutrino events are stable over energy 
and time, and the inferred external event rates are similar 
in both neutrino and antineutrino modes. Furthermore, 
any single background would have to be increased by 
more than 3 a to explain the observed excess of events. 
An additional check comes from the data in neutrino 
mode, which has a similar background to antineutrino 
mode and where good agreement is obtained between the 
data and Monte Carlo simulation for E® E > 475 MeV. 
As a final check, the event rate of candidate v e events in 
the last 2.27 x 10 20 POT is found to be 1.9a higher than 
the candidate event rate in the first 3.39 x 10 20 POT Q; 




10 2 



10 



> 

CD 



E 

< 



10" 



10" ; 



90% CL unsubtracted 
99% CL unsubtracted 
90% CL subtracted 
99% CL subtracted 
KARMEN2 90% CL 
BUGEY 90% CL 




10"' 



sin 2 (26) 



10* 

sin 2 (26) 



10" 



FIG. 3: MiniBooNE 68%, 90%, and 99% CL. allowed regions 
for events with E„ > 475 MeV within a two neutrino v^ — > 
v £ oscillation model. Also shown are limits from KARMEN 
[la] and Bugey [25|]. The Bugey curve is a 1-sided limit for 
sin 2 26» corresponding to Ax 2 = 1.64, while the KARMEN 
curve is a "unified approach" 2D contour. The shaded areas 
show the 90% and 99% CL. LSND allowed regions. The black 
dot shows the best fit point. 

however, the v^ event rates are found to be similar for 
the two running periods. 

Fig. Q] (bottom) shows the event excess as a function of 
E~} E . Using a likelihood-ratio technique, the best Mini- 
BooNE oscillation fit for 475 < E$ E < 3000 MeV occurs 
at {Am 2 , sin 2 20) = (0.064 eV 2 , 0.96). The energy range 
E® E > 475 MeV has been chosen for the fit as this is the 
energy range MiniBooNE used for searching for neutrino 
oscillations. Also, this energy range avoids the region of 
the unexplained low-energy excess in neutrino mode |3J. 
This best-fit point has only a slightly lower x 2 than other 
points in the allowed band. The % 2 for the best-fit point 
in the energy range of 475 < E® E < 1250 MeV is 8.0 for 
4 DF, corresponding to a x 2 -P r °bability of 8.7%. The 
probability of the background-only fit relative to the best 
oscillation fit is 0.6%. Fig.[3]shows the MiniBooNE 68%, 
90%, and 99% CL. closed contours for v^ — > v e oscilla- 
tions in the 475 < E® E < 3000 MeV energy range, where 
frequcntist studies were performed to determine the CL. 
regions. The allowed regions are in agreement with the 



FIG. 4: MiniBooNE 90% and 99% CL. allowed regions for 
events with E„ > 200 MeV within a two neutrino £/ M — > v e 
oscillation model. The solid (dashed) curves are without 
(with) the subtraction of the expected 12 event excess in the 
200 < E„ E < 475 MeV low-energy region from the neutrino 
component of the beam. Also shown are limits from KAR- 
MEN 4j] and Bugey [11 . The shaded areas show the 90% 
and 99% CL. LSND allowed regions. The black dots show 
the best fit points. 

LSND allowed region. The MiniBooNE closed contours 
for Vf, ->• v e oscillations in the 200 < E$ E < 3000 MeV 
energy range are similar, as shown in Fig. 01 The solid 
(dashed) curves are without (with) the subtraction of the 
expected 12 event excess in the 200 < E® E < 475 MeV 
low-energy region from the neutrino component of the 
beam. The best oscillation fits without and with this sub- 
traction occur at (Am 2 , sin 2 20) = (4.42 eV 2 , 0.0066) and 
(4.42 eV 2 , 0.0061), respectively, while the corresponding 
^-probabilities in the 200 < E$ E < 1250 MeV energy 
range are 10.9% and 7.5%. 

A further comparison between the MiniBooNE and 
LSND antineutrino data sets is given in Fig. [5J which 
shows the oscillation probability as a function of L/E u 
for i/fj, — >• v e candidate events in the L/ E v range where 
MiniBooNE and LSND overlap. The data used for LSND 
and MiniBooNE correspond to 20 < E v < 60 MeV and 
200 < E$ E < 3000 MeV, respectively. The oscillation 
probability is defined as the event excess divided by the 
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FIG. 5: The oscillation probability as a function of L/E„ E for 
Vp — > v e candidate events from MiniBooNE and LSND. The 
data points include both statistical and systematic errors. 

number of events expected for 100% v^ — > v e transmu- 
tation, while L is the reconstructed distance travelled 
by the antincutrino from the mean neutrino production 
point to the interaction vertex and E v is the recon- 
structed antineutrino energy. The Lj ' E v distributions for 
the two data sets are consistent. 

In summary, the MiniBooNE experiment observes an 
excess of v e events in the energy region above E® E of 475 
MeV for a data sample corresponding to 5.66 x 10 20 POT. 
A model-independent hypothesis test gives a probability 
of 0.5% for the data to be consistent with the expected 
backgrounds in the energy range of 475 < E® E < 1250 
MeV, and a likelihood-ratio fit gives a 0.6% probabil- 
ity for background-only relative to the best oscillation 
fit. The allowed regions from the fit, shown in Fig. [31 
are consistent with v^ — > v e oscillations in the 0.1 to 1 
eV 2 Am 2 range and consistent with the allowed region 
reported by the LSND experiment [l|. Additional run- 
ning in antineutrino mode is expected to approximately 
double the current number of POT. 
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